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ABSTRACT 



Direct application of two-phase heat transfer in the liquid cooling of electronic 
components in llouorinated hydrocarbons (FC-72), is severely inhibited by the excessive 
amount of superheat required to initiate nucleate boiling. To conduct an experimental 
study of nucleate pool boiling of FC-72, an experimental test chamber was constructed. 
This chamber utilized five horizontal platinum wires of 0.05 mm diameter spaced 2.0 
cm vertically from each other. The lowest wire was progressively heated from the na- 
tural convection region through nucleate boiling, and a study was made on the effects 
of the boiling wake plume on the heat transfer rate of the upper wires. 
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I. INTRODUCTION 



A. HISTORICAL 

Ever since the development of the first practical solid state integrated circuit in 1959, 
there has been a continuous emphasis on packaging an ever higher density of logic cir- 
cuits on a single substrate. Corresponding with this increase packaging density is an 
increasing power dissipation requirement along with an increase difficulty in the thermal 
management capabilities. Initially, thermal management techniques consisted of noth- 
ing more than a heat sink and natural convection air cooling. This naturally progressed 
into heat sinks with forced convection air cooling, and with the construction of large 
mainframe computers and first generation super-computers the use of indirect liquid 
cooling was initiated. The next practical method for heat removal is direct liquid cooling 
or immersion cooling of microelectronic components using dielectric liquids. 

Although direct liquid cooling has been successfully employed for nearly 40 years in 
large hi-powered microwave equipment, klystron tubes, travelling-wave guide tubes, and 
power transistors [Ref. 1 : p. 1 2], only recently has it been employed for use in the ther- 
mal management of microelectronics. The first computer systems to successfully employ 
immersion cooling of microelectronics were the CRAY-2 and the ETA- 10 supercom- 
puter [Ref. 2 : p. 8]. Today the research in direct immersion cooling is continuing with 
particular emphasis on utilizing dielectric liquids as a cooling medium, in particular the 
preflu orinated FC liquid family of the 3M corporation [Ref. 3]. 

B. BACKGROUND 

The advantages offered by direct immersion cooling is self evident from Figure 1 
By utilizing forced convection direct liquid cooling with FC-77, the CRAY-2 has been 
able to reduce it's required volume by an order of magnitude as compared to it's indirect 
liquid cooled competitors. This was being achieved in the CRAY-2 while at the same 
time maintaining the same power dissipation of the competition, [Ref. 1]. 
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Figure 1. Power consumption vs Volume: Of present computers 

As previously stated, the CRAY-2's packaging density was achievable due to direct 
liquid immersion cooling, however, it did use a forced convection cooling scheme. A 
thermal management scheme that relies on boiling heat transfer would provide an even 
more efficient heat transfer capability, and has been the attention of considerable re- 
search in recent years. Unfortunately, these prefluorinated dielectric fluids possesses 
anomalies which have held back their widespread use and acceptance of direct immersion 
cooling. 

In comparison to water and many other coolants, these prefluorinated dielectric 
liquids posses a lower value for their thermal conductivity, and heat of vaporization, 
they also have an extremely small surface tension which translates into a near zero wet- 
ting angle and therefore, a high wettability on most surfaces. This last feature results in 
the elimination of many potential nucleation sites which in turn results in a incipience 
super heat excursion which results in a very high wall temperature. 

In this study, heat transfer enhancement, and therefore a surface temperature re- 
duction, due to bubble pumping in FC-72 will be investigated. 

C. PREVIOUS WORK 

An extensive amount of research has been conducted in the boiling regime of 
dielectric liquids. The most helpful of the literature used will be noted. 
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You, Simon and Bar-Cohen [Ref. 4 and Ref. 5 J conducted an analyses on the ef- 
fects of pressure, subcooling, and dissolved gas content on the boiling heat transfer, and 
the boiling incipience of dielectric fluids. In their experiment they used FC-72 as the 
working fluid and a 0. l^m thin film of platinum on a 0.51mm diameter quartz cylinder, 
with the platinum film acting as both heater and temperature sensor. They concluded 
that: 

1. An increase in pressure of 56kPa reduces wall superheat by approximately 3°C in 
the fully developed nucleate boiling region. 

2. Subcooling has little effect on the fully developed nucleate boiling region. 

3. Nucleate boiling hysteresis is observed for both the subcooled and gas-saturated 
cases. The hysteresis decreases as the gas content increases. 

4. Incipience bubble radius computed for FC-72 was 0.05-0. 1/rm. 

5. Effects of dissolved gases in the fluid is seen at low content values as an increase 
in run to run variability. 

6. The effect of subcooling on pool boiling incipience is apparently small. 

7. A wide variation in incipience superheat values for nominally identical runs were 
observed. 

You, Simon, Bar-Cohen and Tong [Ref. 6] conducted an experimental investigation 
of boiling in R- 113. The experiment was conducted using 0.13mm diameter Chromel 
wires and 0.51mm diameter platinum thin film heater. The following was found. 

1. A Scanning Electron Microscope (SEM) analysis was conducted on the surface of 
the thin film heater, in which the chromel wire surface was found to contain a 
higher density of nucleation sites compared to the thin film platinum heater. 

2. The platinum heater surface was found to have a higher incipience superheat than 
the chromel wires , 73° C versus 38° C maximum. This tends to correlate well with 
item 1. 

3. Again, a wide variation in incipience superheat values for nominally identical runs 
were experience. 

Marto and Lepere [Ref. 7] conducted experiments with surface enhancements for 
improving pool boiling performance. Marto and Lepere using FC-72 and R-113 as 
working fluids compared the performance of three commercially available boiling en- 
hancement surfaces to plain copper tubing. In each case the enhanced surface resulted 
in a lower wall surface temperature compared to the plain copper tube. However, the 
required incipience superheat for nucleation remained unaltered. 

Bergles and Kim [Ref. 2] analyzed methods to reduce the temperature overshoot of 
immersion cooling. Their method used an active bubble generator in the fluid located 
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below the simulated micro chip. They used the bubble generator to activate dormant 
and extinct nucleation sites. They obtained a reduction of incipience superheat required 
from 22-33° C to 8-12° C in a working fluid of R- 113. 

Maddox and Mudawar [Ref. 8] studied the effect of subcooling and surface aug- 
mentation on values of Critical Heat Flux (CHF). They found that the CHF of FC-72 
can achieve a level greater then 360 Watts/cm 2 with combined use of subcooling and 
surface augmentation. 

Danielson, Tousignat and Bar-Cohen [Ref. 9] conducted an analysis of a large 
number of commercially available prefluorinated inert fluids. The method of investi- 
gation used a 0.25mm diameter platinum wire heating element, with the working fluid 
at saturation conditions. Important points reported are: 

1. The critical heat flux of FC-72 was reported as approximately 20.3 Watts/cm 2 . 

2. The nucleate boiling characteristic of highly wetting fluids can be predicted using 
correlations developed for water and other fluids, but are highly dependent of the 
heating surface condition. 

3. Single phase convection data are predicted accurately by a correlation developed 
by Kuehn and Goldstein.[Ref. 10] 

4. The superheat excursion at incipience boiling always appears, but the magnitude 
varies widely. 

5. The apparent radii of nucleation bubbles were calculated to be 0.1 to 0.2/rrn. 

Kuehn and Goldstein [Ref. 10] has obtained a correlation of natural convective heat 
transfer for infinite cylinders. This correlation was used in this thesis data analysis due 
to the low rayleigh numbers obtained in the experimental runs. 

From the 3M Product Manual [Ref. 3] it can be found that the highly wetting 
dielectric fluids have a high solubility for gases. FC-72 can contain up to 48% air by 
volume compared to 1.9% air for water. This is the necessity of thoughly degassing the 
fluid prior to the experimental runs. 

D. OBJECTIVES 

The objective of this work was to investigate the effects of the boiling wake from a 
heated wire upon the heat transfer rate from wires placed within the wake. This work 
was carried out using the dielectric fluid FC-72. 
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II. EXPERIMENTAL APPARATUS 



A. DESCRIPTION OF COMPONENTS 

The experimental apparatus used in this study is a redesign of an experimental ap- 
paratus used in a previous thesis. Emphasize on this design was to minimize systematic 
errors, minimize uncertainty, and to correct deficiencies to ensure an overall increase in 
the data's creditability over the previous apparatus, [Ref. 11]. 

In Figure 2 an overall view of the experimental setup is shown, identifying all major 
components that were utilized. Figure 3 is a photograph of the actual experimental 
chamber. 

The apparatus chamber consists of a box constructed of 0.75 inch thick 
polycarbonate sheets with a 0.25 inch thick aluminum coverplate fitted with a rubber 
O-Ring seal. Inside the box are two bulk heaters to raise the fluid to saturation tem- 
perature for degassing purposes, and to maintain the bulk fluid temperature near con- 
stant. The aluminum cover plate is mounted with three thermo-electric coolers, and 
serves as a condenser, where by the vapor condenses and returns to the chamber by 
gravity. An insert board is held in the chamber on which are mounted the test wires, 
tension springs, thermocouples, and voltage pickofF s. 

1. Test Chamber 

The interior dimensions of the chamber are 6.0 inches wide, by 2.375 deep, by 
6.25 inches high , as shown in the isometric view of Figure 8. The entire chamber was 
assembled using IPN #16 Acrylic bonding adhesive and a bead of Di-Metheyl Ethelene 
was used to seal all seams along the interior surface to ensure liquid tightness. 

To permit access for instrumentation and power leads, and as a means to ensure 
pressure equalization, a 0.75 inch diameter and 0.25 inch diameter hole were placed into 
the rear face of the chamber at the lower right and upper right comer respectively, as 
shown in Figure 4. Mounted securely around these access holes are lengths of tygon 
tubing of approximately 18 inches in length, with the opened ends suspended well above 
the chamber. 

On the interior side walls, slots are milled to permit the placement and removal 
of an instrumented wire board on which the experiments are conducted. Figure 6 is a 
photograph of the insert board, while a schematic drawing is shown in Figure 7. 
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2. Aluminum Cover Plate 

A 0.25 inch thick by 7.00 inch long by 3.25 inch wide aluminum plate was uti- 
lized for use as a cover plate for the chamber, and as a condenser surface for the FC-72 
vapor. Mounted on the interior face are two thermocouples to monitor the condensor 
surface temperature, as shown in Figure 9. On the exterior face of the plate, three 
thermo electric coolers were installed for heat removal. To ensure a proper seal between 
the aluminum plate and the box structure , a shallow O-Ring groove was machined into 
the interior face and an O-Ring of 0.0624 inch diameter was installed. To adequately 
secure the cover plate to the chamber box, eight stainless steel screws of 1/16 inch di- 
ameter are used to compress the O-Ring against its mating surface. 

As mentioned previously, three solid state thermoelectric cooling devices from 
the MELCOR corp were employed to cool the cover plate condensor. They were nom- 
inally operated about the 0.4 amps and 2.0 volts operating point, from [Ref 12] their 
operation curves and technical data is available. 

3. Heaters 

As a means of degassing the FC-72 liquid and for bulk temperature control, two 
strip heaters of 125mm x 1 1 mm were employed. These strip heaters were secured to the 
bottom interior surface of the chamber with epoxy, and were orientated lengthwise in the 
chamber such that one heater strip will be in front of and one heater strip behind the 
insert board. To raise the bulk fluid temperature to saturation temperature and for 
degassing purposed, both heaters are on line and operated at 1.8 amps and 20 volts. 
Once degassing is complete, the forward strip heater is electrically isolated and the rear 
heater is run at 0.9 amps and 20 volts. This will maintain the forward portion of ex- 
perimental chamber at near saturation temperature with minimum stratification of fluid, 
and maintain liquid in a degassed state since the fluid is still boiling from the heater 
strip's surface. 

4. Thermocouples 

A total of Six copper-constantan type thermocouples with a wire diameter of 
0.05 inches (0.127 mm) were utilized for temperature measurements. Two thermocouples 
were bonded to the underside of the aluminum condensor plate utilizing OMEGA Bond 
type 101 thermally conductive epoxy, four thermocouples were installed on the insert 
board, and one was used for measuring the ambient air temperature. 

5. Insert Board 

An insert board was designed to mount into two 0.250 inch wide grooves that 
were milled into the chamber side walls. Onto this insert board are mounted four 



6 



thermocouples to measure bulk fluid temperature, five platinum wires of nominally 0.002 
inch (0.05mm) diameter by 4.0 inch long for conducting the actual boiling experiments. 
Five springs are used to maintain the platinum wire tension so no sagging develops due 
to thermal expansion of the wires. In Appendix A the platinum wire calibrations and 
specifications are shown in detail. 

B. INSTRUMENTATION 

In the design and assembly of this experimental apparatus, five channels were wired 
and instrumentation provided to measure voltage drops at select points. Each of these 
channels had a 2.0 ohm precision resistor wired in series with the platinum wire, where 
the voltage drop across the resistor was used to obtain the current flow through the 
platinum wire. The voltage drop measured across the platinum wire is used to obtain 
the wire resistance, power dissipated and the heat flux of the platinum wire. The 
platinum wire temperature is obtained from its resistance vs temperature calibration 
curve obtained prior to experimental runs. Instrumentation was additionally used to 
obtain the bulk fluid temperature, which will be required for calculating fluid properties. 

1. Platinum Wire 

Each of the five platinum wires were run in series with a 2 ohm (± 1%) precision 
resistor of 25 watt power rating. The following table gives the platinum wire and resistor 
combinations which were used. 



Table 1. MEASURED RESISTANCE VALUES FOR 2 OHM RESISTORS 



Wire Chan- 
nel 


Measured Resistance, 20 Sample 
average 


20 Sample Std. Dev. 


1 


1.9994^ 


0.0009a 


2 


2.0076a 


0.0005a 


3 


1.9998a 


0.0003a 


4 


2.0018a 


0.0013a 


5 


2.0025a 


o.ooosa 



A 2.0 ohm precision resistor was chosen over the previous design's 0.1 ohms 
precision resistor in order to reduce the calculated currents uncertainty, and there by 
reduce the calculated platinum wire temperature uncertainty, [Ref. 11: p. 25]. 
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2. Power Supplies 

The following power supplies were utilized and their primary function listed. 

1. Hewlett Packard 6214C 0-10volts/0-l amps: Used to power the thermal electric 
coolers. 

2. Hewlett Packard 6214C l-10volts / 0-1 amps: Used to power the constant heat flux 
platinum wire. 

3. WP71 1 0-40 volts / 0-1 amps : Used to power the varying heat flux platinum wire. 

4. Hewlett Packard 6286A 0-20volts / 0-10 amps : Used to power Heater Strips for 
degassing and bulk temperature maintenance. 

5. Hewlett Packard 6289A 0-40 volts / 0-1.5 amps : Used to power an auxiliary illu- 
mination system or an augmenting cooling fan for the thermal electric coolers. 

3. Acquisition Unit 

To obtain the necessary measurements, a Hewlett Packard 3852A Data Acqui- 
sition System was employed, utilizing a Hewlett Packard 300 series desktop computer. 
The HP3852A utilized the following utility modules inserted into it's chassis. 

1. Hewlett Packard Model 44701A Digital Integrating Voltmeter 

2. Hewlett Packard Model 44705A 20 channel relay multiplexer for measuring re- 
quired voltage drops. 

3. Hewlett Packard Model 447 13A 24 channel relay multiplexer with electronic cold 
junction for direct conversion of thermal couple emf to temperature. 

The HP300 desktop computer was used to control the HP3852A Data acquisition sys- 
tem. Additionally, through specifically written programs, the computer can accomplish 
the required data reductions, (See Appendix D). 



8 



h 



X. 




r 



\ 

\ 

\ \ 



o lli 
3 I- 
d LO 
o >- 

< (O 



cr 

LU 

h- 

3 

CL 

o 

o 



OvJ 

I 

CO 

CL 



x\ 


o_ 


s' 


— 


/ 


o 



(O 

Q_ 



\ 



/ 



o. 

o' 



! I 



cr 

LU 
1 -CD 
(02 
LU < 

>-x 
\ o 
\ 



\ 



“1 



1 






\N s 



N 



\ 




L(~)| fO C\J — ^ 






V 



X 



r 



m 



m 

I 

LD 

CL 



vs 

5 



o_ 

o~ 



_ i 



Figure 2. Overall experimental settup 
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Figure 3. Experimental Chamber in Detail 
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Figure 4. Side View of Experimental Apparatus 
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Figure 5. Rear View of Experimental Chamber 
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Figure 6. Insert Board 
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Figure 7. schematic Drawing of Insert Board 
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Figure 9. Condensor Surface Details 
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III. EXPERIMENTAL PROCEDURE 



A. PREPARATIONS 

Careful preparations were required to ensure that the data obtained was of high 
quality. The most notable requirement being the necessity of thoughly degassing the 
FC-72 liquid, [Ref, 3: p. 3]. 

1. Normal Warm-Up Procedure 

In carrying out the experimental runs the following procedures were followed. 

1. Power-up the HP3852A Data Acquisition Unit to start the warm-up period of one 
hour minimum. [Ref. 13: p. 586] 

2. Ensure the liquid level is above Wire 5, but below the condensate drip tray. 

3. Set power supply of thermoelectric coolers to 1.3 amps and 2.0 volts. 

4. Ensure the experimental chamber's insulation jacket is in place and set power sup- 
ply to strip heaters to approximately 1.8 amps and 20 volts. 

5. Verify the 0.25 inch diameter vent line is unobstructed. 

6. Power on the HP300 series desktop computer, and load data acquisition program 
as contained in Appendix D. 

7. Once fluid has reached saturation and is boiling vigourisly, hold for a minimum of 
one hour, typically two hours was achieved, to permit adequate degassing of fluid. 

In actual practice, the time required to raise the bulk temperature of the fluid to satu- 
ration temperature of 56 °C was nominally two hours. 

2. Post Warm-up Procedure 

Once the required degassing has been accomplished, power is temporarily se- 
cured to the strip heaters. Once power is secured the forward strip heater is electrically 
isolated. The rear strip heater is powered at 0.9 amps and 20 volts, and the thermo 
electric coolers are secured. 

The power setting of the strip heaters permits vigorous boiling from it's surface. 
This serves to keep the fluid degassed and the bulk fluid temperature near saturation and 
near isothermal. By keeping the fluid level below the insert board's upper edge, the fluid 
forward of the insert board remains undisturbed by the boiling action occurring in the 
rear chamber. 
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B. DATA ACQUISITION PROCEDURE 

Upon completion of the previous preparatory procedures, the experimental runs are 
ready to begin. The procedure followed for an experimental run is as follows: 

1. Plug power supply into power receptacles for wire 1. Power on and ensure current 
limiter is set to max value of 1.0 amps. 

2. Plug HP6214C power supply into power receptacles of desired constant heat flux 
wire. 

3. Set power supply of step 2 to proper output to give desired heat flux from the wire. 

4. Start running the Data Acquisition Program, which is programmed to measure the 
following data from the desired wires at two minute intervals. 

• Precision Resistor Voltage drop, up to five channels. 

• Platinum Wire Voltage drop, up to five channels. 

• Four bulk fluid temperature measurements 

• Two condenser surface temperature measurements. 

• One ambient air temperature measurement. 

5. After each data interval, increase the voltage to the varying heat flux wire by 0.2 
Volt increments. This was continued until the on line power supply was at 2.3 volts 
of outputs, the increments were then reduced to 0.1 volts until nucleation occurred. 
The increments were then returned to 0.2 volts until a heat flux of approximately 
200,000 watts/m 2 were achieved. At any higher heat flux the wire will go into the 
film boiling regime. 

6. Reduce the power to the varying heat flux wire at 0.3 Volts increment to bring the 
wire back to a no power condition. 

7. Repeat steps 3 thru 6 for all other wire and heat flux configurations. 



C. DATA REDUCTION 



Data that was obtained from the actual experimental apparatus were two voltage 
measurements per active platinum wire channel, and six temperature measurements. 
The process used for data reduction is a combination of computer programs and 
graphical outputs for numerical presentation and visual presentations respectively. 

1. Numerical Data 

The numerical data values required to be processed are: 

1. Fluid Bulk Temperature 



(T,+T 2 + T 3 + T 4 ) /o _ 

1 bulk — l M 



(3.1) 



and where, T, thru T 4 are the thermocouple measurements of the bulk fluid. 
2. Current flow through the Platinum Wire, where 
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1 wire = 4 5SL (Amps) 
K 2 Q 



3. Platinum Wire's electrical resistance value, where 

V PtWlre /r\\- \ 

^ Pi Wire i (Ohms) 

1 wire 

4. Platinum Wire surface temperature, where 

Tmr = Z^Pimre + T 0 (°C) 

where, 

T,„, = the actual wires surface temperature 
e is the slope of the calibration curve 
T„ is the calibration curves y-axis intercept 

5. Fluid Film Temperature, where 

T _ + Jbulk 

Iy?/m - 2.0 



6. Dissipated Heat Flux of platinum wire, where 



q = 



1 wire ^ y PtlVire 



Area 



Surf 



, Watts .. 

( 9 ) 



m 



where, Area w = n x length.,,,, x D w/ „ 



(3.2) 



(3.3) 

(3.4) 



(3.5) 



(3.6) 



The above numerical data were all derived from the raw data, and previously 
determined calibration curves. The following table provides physical dimensions and 
temperature calibration coefficients of the platinum wires. 



Table 2. PLATINUM WIRE PHYSICAL PARAMETERS FOR EQUATION 



(3.4) AND EQUATION (3.5) 



Wire 


6 


T, 


length 


diameter 




o 

n 

B 


°C 


mm 


mm 


1 


53.902 


-254.470 


99.22 


0.05 


2 


53.666 


-253.009 


99.22 


0.05 


3 


53.141 


-252.729 


100.01 


0.05 


4 


53.484 


-252.848 


98.43 


0.05 


5 


53.086 


-260.049 


100.01 


0.05 



To see in detail the methodology used to obtain the numerical values presented 
in the preceding table, proceed to Appendix A 'Platinum Wire Calibration'. 
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2. Fluid Property Data 

The Fluid Properties to be determined are as follows: 

1. Thermal Conductivity of fluid, as obtained from Figure 5 of the 3M Corporation 
Flourinert Product Manual [Ref. 3: p. 15] have been determined to be, 



0.6033 — 0.001 15 x Ty^ Watts x 

k la 



(3.7) 



2. Thermal Expansion Coefficient, from table 4B in the Product Manual [Ref. 3: p. 10] 
is given as 



„ 0.00261 , 1 . 

P 1.740 — 0.00251 x Tyj /m °C 



(3.8) 



3. Kinematic Viscosity of the Fluid , as obtained from Figure (3) in the Product Man- 
ual [Ref 3: p. 1 3] and determining an exponential curve fit yields 



v= 1.203952E-8 x exp(l, 058.4109° 



K 



T, u/ * + 273.15°C 



X^) 



(3.9) 



4. Specific Heat of the Fluid for all the Flourinert Fluids are obtained from Figure (4) 
of the Product Manual [Ref 3 : p. 14] is, 



C, = (0.24111 + 3.70374£— 4 x T bulk ) x 4186( 



kg°C 



) 



(3.10) 



5. Liquid Density of FC-72, from expression in table 4B and constants from table 4C 
in the Product Manual [Ref 3: p 10] is given by 



p = (1.740 - 0.00261 x T to/fc ) x 1000( ) 

m 



6. Thermal Diffusivity of the liquid is given by, 

k 



a = 



pxC. 



(3.11) 



(3.12) 



7. Prandlt Number of the Fluid is given by 

Pr = -F 

8. Ralyeigh Number of Fluid based on wire diameter is given by, 



Ra = 



g x fi x (T wr - T x D 
v x a 



(3.13) 



(3.14) 



9. Correlation for the Nusselt Number was obtained by using the correlation of Kuehn 
and Goldstein [Ref 10: p. 1128] given by 



Nuq = 



ln[l +■ 



[(0.518 x Ra if x [1+ ( yi>y 3 i i y 3 + (o.l x Ra)j J ) ,P ] 



0.559 ^3/5-,— 5/12. 15 



1/3.15^1/15 
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10. Theoretical Heat Flux as obtained from using the N D value obtained from equation 
(3.15). 



k x Nu D x (T rar - T bu!k ) Watts . 

' 2 ' 



D 



m 



(3.16) 



With all the required numerical data and fluid properties obtained, it is now 
necessary to present the numerical information in a graphical representation to see the 
relevancy of the results obtained. This will become apparent with the presentation of 
results in Chapter IV 'Results and Discussion'. 
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IV. RESULTS AND DISCUSSION 



The data obtained can be grouped for discussion as follows: 

1. Individually powered wires 

Saturated Conditions 
Sub-Cooled Conditions 

2. Bubble Pumping Effect Runs 

Saturated Conditions 
Sub-Cooled Conditions 

A. INDIVIDUALLY POWERED WIRES 

Before analyzing the pumping effect of the bubbles on the upper wires, the individ- 
ual wire boiling experiments will be investigated. All of the individual powered wire runs 
shown here were made at or near saturated fluid conditions with two exceptions. The 
following table will provide a listing of single wire runs. 



Table 3. INDIVIDUAL WIRE RUNS 



Figure 


Wire 


T bulk 


Figure 10 


1 


55.5° C 


Figure 11 


1 - several runs overlayed 


At or near saturation 


Figure 12 


1 


38° C 


Figure 13 


1 


23.0° C 


Figure 14 


2 


56.0° C 


Figure 15 


3 


56.0° C 


Figure 16 


4 


56.0° C 


Figure 17 


all the above 


see individual wires 



In all cases conditions were under atmospheric pressure with 0.05mm diameter platinum 
wires. 

1. Natural Convection 

Figure 10 is a typical boiling curve for FC-72 with the Kuehn and Goldstein 
natural convection region shown along with the experimental points. The nucleate 
boiling zone starts at boiling incipience, which is evident from the sudden drop of wall 
temperature, a characteristic of a highly wetting fluid. Figure 1 l'is a compilation of wire 
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1 data, at or near saturation, obtained over a period of time. As can be seen, with the 
exception of the boiling incipience point, the data is virtually super imposable over each 
other. This reproducability of data over time indicates no significant aging has occured 
on wire 1, and more importantly no change in calibration. 

In fact in the natural convection region Figure 10 thru Figure 16 are all essen- 
tially superimposable over each other. Figure 17 is an overlay of all the previous fig- 
ures. Therefore, it appears that whether the fluid is Saturated or Sub-Cooled, the fluid 
bulk temperature appears to have little effect on the natural convection heat transfer 
rate. 

There was generally very good agreement between the natural convection data 
and the Kuehn and Goldstein correlation [Ref 10] (equation 3.15). The only deviation 
existing between the figures being the total temperature excursion to incipience, which 
ranged from 27° C to 38° C. 

2. Nucleate Boiling 

As previously stated, the nucleate boiling zone starts at boiling incipience. Us- 
ing wire 1 as a standard and Figure lOas base line boiling characteristic data a compar- 
ison to other wires at saturation condition, and a comparison to subcooling effects on 
wire's 1 boiling can be made. As is evident from Figure 17, the boiling characteristic 
of all the wires at saturation condition tend to overlap each other. The differences aris- 
ing tending to be caused by individual wire placement in the fluid . However, the effect 
of subcooling does result in a decrease in the boiling curve overshoot at point of 
incipience. This decrease in the boiling curve overshoot results in a net decrease in the 
superheat required to initiate nucleate boiling. 

3. Departure from Nucleate Boiling 

The Critical Heat Flux (CHF) for FC-72 at saturated condition was calculated 
by Danielson et all [Ref 9] to be approximately 203,000 Watts/m 2 . Experimental value 
for CHF ranged from 198,000 Watts/m 2 to 220,000 Watts/m 2 , a spread of -2.5% to 
8.4%. A CHF for a subcooled condition of T Aa/A =38° C was found experimentally to 
be approximately 400,000 watts/m 2 . A CHF for T 4u/A = 23 0 C was not obtained during 
the experimental runs. 
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WIRE 1 BOIUNG CURVE 



NEAR SATURATION T bulk = 55.5 °Celisus 




Figure 10. Wire 1 boiling curve T 4u/ * = 55.5°C 
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WIRE 1 BOILING CURVES 



ALL BOILING CURVES AT SATURATION OR NEAR SATURATION 




Figure 11. 



Wire 1 Boiling Curves over time 
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WIRE 1 BOILING CURVE 




Figure 12. Wire 1 boiling curve T iulk = 38°C 
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WIRE 1 BOILING CURVE 



.SUBCOOLED T buik = 22.5° Celisus 
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Figure 13. Wire 1 Boiling Curve T 6u/ * = 23°C 
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WIRE 2 BOILING CURVE 



SATURATED T bu , k = 56.0°Celisus 




Figure 14. Wire 2 boiling curve T w * = 56°C 
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WIRE 3 BOILING CURVE 



SATURATED T Mk = 56.0° Celisus 




TSURF - TSAT (C) 



Figure 15. Wire 3 Boiling Curve T 6urt = 56°C 
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WIRE 4 BOIUNG CURVE 



SATURATED T bulk = 56.0° Celisus 




TSURF - TSAT (C) 



Figure 16. Wire 4 boiling curve T tu/ * = 56°C 
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BOILING CURVES COMPILED 
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Figure 17. Overlay of previous curves 
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B. BUBBLE PUMPING EFFECT ON HEAT TRANSFER RATE 



In the upcoming discussion a brief explanation of the figures are required. In each 
case for these experiments, heater wire number 1 was used to generate the bubbles for 
examining the pumping effects upon the upper wires. The figure insert is used to show 
which wires are being utilized. Only select data points were plotted for the passive 
(constant heat flux wire) wire, and those points correspond in time to the identically 
numbered data point on wire l's boiling curve. The data points were chosen since they 
define regions of wire l's boiling curve as shown in the following table. 



Table 4. DATA POINT DESCRIPTION 



5 point figure 


6 point figure 


Region Description 


1 


1 


Natural Convection Start 


2 


2 


Point of Incipience 




3 


Partial Boiling on wire 


3 


4 


Full nucleate Boiling 


4 


5 


Max Nucleate Boiling 


5 


6 


End of run 



1. Passive Wire at Low Heat Flux 

In these experimental runs the heat flux on the passive wire was set to approx- 
imately 8500 watts/m 2 . At this low heat flux it was expected that little to no active 
nucleation will occur on the passive wire, even after wire 1 is in nucleate boiling. 
Therefore, all the enhancement of the upper wire will be due to an addition of a force 
convection term to the heat transfer rate. 

In comparing Figure 18 to Figure 19the following can be inferred: 

1. Total reduction of wall temperature on the passive wires are comparable. Wire 2 
having a 3.7°C and Wire 4 having a 3.0° C drop in wall temperature. 

2. A definite effective distance exist over which wire's 1 influence is imparted upon the 
upper wires, for example: 

Wire 2 and wire 4 showed a very small change in wall temperature at point 2 (point 
of incipience). 

After incipient boiling has occurred on wire 1, wire 2 showed a greater drop in wall 
surface temperature then wire 4. Approx. 2.5° C drop for wire 2 and 0.1° C drop 
for wire 4. 
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3. At point 4, the point of maximum temperature in nucleate boiling, both wire 2 and 
wire 4 showed a comparable surface temperature, approximately 2° C above satu- 
ration. 

As can be seen from the figures, some nucleation sites have been activated on wire 2 and 
wire 4, with point 5 showing a lower wall temperature then point 1. 

2. Passive Wire at Medium Heat Flux 

In these experimental runs the heat flux on the passive wire was set to approx- 
imately 20,500 watts/m 2 . At this heat flux it was expected that active nucleation will 
occur on the passive wire after wire 1 is nucleate in boiling. 

In comparing Figure 20 to Figure 21 
the following can be inferred: 

1. The total reduction of wall temperature on the passive wire , points 1 to points 4, 
are comparable at 7.1° C for wire 2, and 6.3° C for wire 4. 

2. The spatial effect is evident on the influence of wire 1 on the upper wires, for ex- 
ample: 

Between point 1 and point 2, wire 2 shows a greater reduction in wall surface tem- 
perature than is evident in wire 4. 

At point 3, wire 2 shows a much greater drop in surface temperature than wire 4: 
4° C compared to 2° C. The difference being the amount of 

3. At point 4, maximum nucleate boiling from wire 1, both wire 2 and wire 4 show 
comparable values for surface temperature, approximately 

1.3° C. 

As can be seen from the figures, wire 1 has activated previously inactive nucleation sites 
on wire 2 and wire 4. This is evident from point 5 being at a lower surface temperature 
than point 3. 

3. Passive Wire at High Heat Flux 

In these experimental runs the heat flux on the passive wire was set to 40,000 
Watts/m 2 on wire2 and 35,000 Watts/m 2 on Wire 4. At this heat flux it was expected 
that a large number of inactive nucleation sites will be activated after wire 1 is nucleate 
in boiling. 

In comparing Figure 22 and Figure 23 it is evident that the information ob- 
tainable from the figures is comparable to the information for the medium heat flux case. 
However, the wall temperature reductions on the passive wires are much greater as 
compared to the Medium Heat Flux case, in particular the following: 

1. In comparing the surface temperature reduction from point 1 to point 2, wire 2 has 
dropped 2.5° C w’hile wire 4 had almost a zero drop of temperature. This appears 
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to be due to the relative locations of the passive wires with respect to wire 1, and 
is especially apparent at the higher heat fluxes. 

2. Many inactive nucleation sites have been activated by wire l's nucleation bubbles. 
As can be seen by comparing point 1 to point 5(6), a difference on Wire 2 of 13° 
C and on Wire 4 of 7° C is apparent. 

4. Wire 1 Passive with Wire 2 varying 

In these experimental runs, wire 1 was set to 150,000 Watts/m 2 so it was boiling 
vigorously, wire 2 was actively varied over a large heat flux range as can be seen from 
Figure 24. On Figure 24 a previously obtained boiling curve of Wire 2 was overlayed 
onto it. Upon comparing the two curves the most striking feature is the disappearance 
of the Hysterises Loop as is evident in the overlayed boiling curve data. In this ar- 
rangement, the boiling curve overshoot has been removed as a design problem, bubbles 
from Wire 1 has removed a major design consideration. 

5. Subcooled 

In these experiments, the bulk fluid was degassed and then subcooled to a tem- 
perature of 38° C. In these experimental runs, wire 2 was set at 35,000 watts/m 2 and 
53,000 watts/m 2 , as shown in Figure 25 and Figure 26 respectively. In subcooling, no 
significant effects were notable beyond those seen in the saturated runs. The only effect 
that was apparent was the ability to run wire 2 at a greater heat flux with out incipent 
boiling occuring on the wire. 
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HEAT FLUX (WATTS/M2) 



WIRE 2 AT 8,500 WATT/M 2 



NEAR SATURATED T iM = 55.5° Celsisus 




Figure 18. Wire2 at 8500 W/ m J T iutt = 55.5°C 
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WIRE 4 AT 8,500 WATT/M2 



NEAR SATURATION T bu , k = 55.6° Celsisus 




Figure 19. Wire4 at 8500 W/ m J T iurt = 55.6°C 
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WIRE 2 AT 20,500 WATT/M 2 



NEAR SATURATION T iu , k = 55.6°Celisus 




Figure 20. Wire 2 at 20,500 W/ m 2 T w * = 55.9degreeC 
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WIRE 4 AT 20,500 WATT/M2 



NEAR SATURATION T Mk = 55. r Census 




Figure 21. Wire 4 at 20,500 W/ m 2 T w * = 55.6°C 
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WIRE 2 AT 40,000 WATT/M2 



NEAR SATURATION T iulk = 54.6* Celisus 




Figure 22. Wire 2 at 40,000 W/ m 3 T 4urt = 54.6°C 
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WIRE 4 AT 35,000 WATT/M2 



NEAR SATURATION = 55.0° Census 




Figure 23. Wire 4 at 35,000 W/ m J T 4u/ * = 55°C 
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WIRE 1 @ 150,000 WATTS/M 2 



SATURATED T bu , k = 56° C 




Figure 24. Wire 1 passive, Wire 2 active T tutt = 56°C 
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WIRE 2 AT 35,000 WATT/M 2 



SUBCOOLED T bu/k = IS.O'Celsisus 




Figure 25. Wire 2 at 35,000 W / m J T** = 38°C 
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WIRE 2 AT 53,000 WATT/M 2 



SUBCOOLED r 4oft = 38 ° Celisus 




Figure 26. Wire 2 at 53,000 W/ m 2 T iu „ = 38°C 
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V. CONCLUSIONS 



The study conducted in this thesis was to investigate the effects of the wake from a 
bubble generating wire on wires located downstream within the wake. Experimental 
runs were conducted at saturated condition and subcooled conditions in which the fol- 
lowing were found: 

A. SINGLE WIRE BOILING 

1. The effects of subcooling on the natural convective heat transfer from the wire was 
negligible. 

2. The effects of subcooling tended to increase the heat flux required for incipient 
boiling to initiate. 

3. With subcooling, boiling hystersis is still observed, but the hystersis magnitude de- 
creases as subcooling increases. 

4. Subcooling results in a reduction in wall temperature after the start of nucleate 
boiling as compared to saturated cases. 

B. PASSIVE WIRE ANALYSIS 

Heat Transfer enhancement effects from Wire 1 upon the downstream wires has 
many factors within it. The following factors have been identified: 

1. An enhancement to convection resulting from buoyancy induced fluid flow off of 
wire 1 was seen to exist. This effect is very spatially dependent with the effects 
most pronounced on wire 2 and negligible on wire 4. This corresponds to point I 
and point 2 on the figures. 

2. The sudden drop of wall temperature between point 2 and point 3/4 was seen to 
be the results of the following: 

Nucleation bubbles from wire 1 activating the inactive nucleation sites on the upper 
wires resulting in boiling heat transfer off of the upper wires. 

An additional force convection enhancement due to the nucleation bubbles from 
wire 1 breaking up the thermal boundary layer about the upper wires is believed to 
exist. 

3. The drop in temperature on the upper wires corresponding with points 3/4 to point 
4/5 appears to be due purely to an increase in force convection due to the high 
density of bubbles off of wire 1 . 

Additional conclusions reached are. 

1. A constant source of bubble generation from below is an effective way to eliminate 
the Boiling Curve Overshoot and Hysteresis loop associated with dielectric fluids. 
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2. At low heat fluxes, proximity to Wire 1 is important for heat transfer enhancement 
to occur. This dependency reduces as heat flux increase on the passive wire. 
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APPENDIX A. CALIBRATION 



A. PLATINUM WIRE 

1. Surface Micro-Structure 

Prior to the experiments, a length of 0.05mm platinum wire was analyzed using 
a Scanning Electron Microscope (SEM). The length of wire analyzed was not used in 
any experiments, but it did come from the same spool of wire as the experimental lengths 
of wire, therefore, it should be a representative sample. In analyzing this sample the 
SEM brought out several important features: 

1. From Figure 27 a wire diameter of 0.052mm can be determined. Also, the surface 
is relatively unmarked at this magnification of 1490x; little scoring on the wire and 
a low density of surface pits are evident. 

2. In Figure 28, an increase in magnification to 4740x within the field of Figure 27, 
a more detailed view of the larger surface defects are evident. Maximum defect size 
is approx 2^m in length and 0.5^m in width. A surface defect of that size would 
be too large to be a possible nucleation site , according to the calculations done in 
Ref. 4 of 0.1 nm. 

3. In Figure 29, a further increase in magnification to 13,200x within the field of 
Figure 27, pits of approximately 0.8/um max down to 0.3^m are seen. It will be 
these surface defects that boiling will initiate from. The density of these pits along 
the length of the wire analyzed was not great, this implies that the incipient super- 
heat will be large. 






Figure 27. Surface micro geometry: Sample Platinum Wire (1490x) 
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Figure 28. Surface micro geometry: Sample Platinum Wire (4740x) 



Figure 29. Surface micro geometry: Sample Platinum Wire (13,200x) 
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2. Desirability for Using Platinum 

Platinum was chosen for several reasons: 

1 It has a highly linear response of Temperature vs Resistance. 

2. It has a large AR for a given AT. 

In each aspect, platinum performed very well. As can be seen in the calibration curves 
of Figure 31 thru Figure 40 response is linear, with a 0.1£2 change in resistance for ev- 
ery 5° C change in temperature. 

B. BACKGROUND 

An important step in minimizing experimental data uncertainty is an accurate cali- 
bration of the platinum wires used. In obtaining the calibration curves for the wires, 
three calibration runs were performed. 

1. An initial calibration in a heated calibration bath with a platinum resistance 
thermometer. 

2. An in place calibration of wires in the test chamber, with the four installed 
thermocouples used for temperature measurements, was conducted to ensure the 
calibration curves from the calibration bath were not altered from handling the wire 
insert board. 

3. A final in place calibration of wires in the test chamber, 

with the four installed thermocouples used for temperature measurements, was 
conducted to ensure the calibration curves used for data acquisition were not al- 
tered due to any reason. 

The second calibration run was required to ensure no distortion occurred to the wires 
when the insert board, on which the wires are mounted, was removed from the calibrated 
bath and installed into the experimental test chamber. 

C. INITIAL CALIBRATION 

The initial calibration was conducted using the following equipment. 

1. Rosemount Engineering Co. Model 913A calibration bath with Etheylene Glycol 
as the working fluid. 

2. Rosemount Engineering Co. Model 923B power supply, controlling calibration 
bath temperature. 

3. Rosemount Engineering Co. Model 920A commutating bridge, used to measure an 
resistance to 0.001 £2 accuracy. 

4. Rosemount Engineering Co. Model 162C S/N 985 Platinum Temperature Probe. 
Used to measure calibration bath temperature through the commutating bridge. 

5. Hewlett Packard Data Acquisition System. Used to measure platinum wire resist- 
ance using a 4-wire ohm measurement technique. 
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1. Calibration Heating Bath 

The present investigation was conducted using an externally heated calibration 
bath with etheylene glycol as the working fluid. This fluid was circulated by a centrifugal 
pump to ensure a uniform bulk temperature throughout the bath. Item (1) and item (2) 
are utilized . 

2. Temperature Measurements 

Temperature measurements for calibration were conducted using a high sensi- 
tivity commutating bridge to measure the resistance of a platinum resistance temper- 
ature probe. This platinum temperature probe had been previously calibrated to a 
standard that is traceable to the National Bureau of Standards. 

3. Resistance Measurements 

In calibrating the platinum wire the resistance of the wire was measured with the 
same Data Acquisition Unit Components as would be used in the experimental runs. 
This was to ensure the uncertainty of measured resistance would be the same as between 
the calibration run and experimental runs. 

The Data Acquisition Unit Component used are: 

1. HP3852A Control Unit 

2. HP44701A Integrating Voltmeter card 

3. HP44705A 20 channel Multiplexer card 

4. HP300 series desk top computer 

5. HP44713A 24 channel multiplexer with electronic cold junction for thermo couple 
measurements. 

4. Calibration Procedure 

By using the previously described components for measuring temperature and 
resistance, the following procedure was used to obtain the system's calibration data: 

1. Set calibration bath controller for 35° C 

2. Allow bath temperature to stabilize. 

3. Carefully place Insert Board into bath liquid. EXTREME care must be taken to 
ensure that none of the wires are deformed or destroyed. 

4. Allow insert board to come to thermal equilibrium with the bath, 

this takes approximately five minutes. 

5. Take temperature reading with platinum wire temperature probe and commutator 
bridge, and resistance measurement of wires with data acquisition unit. 

6. Remove insert board, and set bath temperature controller for a 5° C increase. 

7. Repeat step 2 thru 6 until an 90° C reading has been achieved. 
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8. Activate External Cooling system. 

9. With temperature control setting decreasing at 5° C intervals, repeat step 2 thru 6 
until a 35° C reading has been achieved. 

D. IN PLACE CALIBRATION 

To ensure that the calibration curve of platinum wire 1 thru 5 are reliable, an in 
place calibration run was performed on the wires. This calibration procedure was per- 
formed before any experimental runs were conducted, and after the last experimental run 
was completed to detect any changes in calibration due to wire deformation or aging. 

The equipment utilized for the inplace calibration was the same Data Acquisition 
equipment used for the resistance measurements of the Bath Calibration runs with two 
additions. 

1. HP44713A 24 Channel Multiplexer with electronic cold junction. 

2. Four type T thermocouples. 

The thermocouples are now utilized, instead of the platinum resistance temperature 
probe, for measuring fluid bulk temperature. The uncertainly associated with the 
thermocouples is 0.4°C as compared to the platinum resistance temperature probe used 
as a standard Figure 30 shows the thermocouples uncertainity compared to the 
platinum resistance temperature probe. 

1. Inplace Calibration Procedure. 

The following procedure is used for inplace calibration of the platinum wires: 

1. Ensure Data Acquisition system is wired for 4-wire ohm measurements. 

2. Turn on Data Acquisition system and allow one hour of warm-up time. 

3. Load Calibration Acquisition Program into the HP-300 series Computer. See Ap- 
pendix D for program listing. 

4. Set power supply to bulk heaters to 20volts and 1.8 amps. 

5. Start Acquisition program. The acquisition program will measure the resistance 
of each platinum wire, and measure the temperature output of each thermocouple 
at a 10 minute interval. 

6. Allow acquisition program to run until bulk fluid has reached Saturation temper- 
ature of 56° C. 

E. CALIBRATION RESULTS 
1. Platinum Wire 1 

From Figure 31 the importance of conducting an in place calibration prior to 
any experimental runs becomes obvious, wire l's Bath calibration data has been con- 
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siderably altered while it was being inserted into the test chamber. The initial in place 
calibration data was used for the experimental runs calibration curve. The resulting 
calibration was: 



T, = 53.902 x R, - 254.470 (°C) 

From Figure 32 it can be seen that the post calibration curve matches well with 
the data acquisition curve. 

2. Platinum Wire 2 

From Figure 33, it is apparent that the initial in place calibration curve fairs in 
well with the calibration bath's decreasing data. The initial in place calibration data was 
used for the experimental runs calibration curve. The resulting calibration was: 

T 2 = 53.666 x R 2 - 253.009 (°C) 

From Figure 34 it can be seen that the Post run calibration matches well with 
the Data Acquisition Calibration curve. 

3. Platinum Wire 3 

From Figure 35 it is apparent that wire 3's Bath calibration data has been al- 
tered between the calibration bath and insertion into the test chamber. The initial in 
place calibration data was used for the experimental runs calibration curve. The result- 
ing calibration was: 



T 3 = 53.141 x R 3 - 252.729 (°C) 

From Figure 36 it can be seen that the Post run calibration matches well with 
the Data Acquisition Calibration curve. 

4. Platinum Wire 4 

From Figure 37, it is apparent that the initial in place calibration curve fairs in 
well with the calibration bath's decreasing data. The initial in place calibration data was 
used for the experimental runs calibration curve. The resulting calibration was: 

T 4 = 53.484 x R 4 - 252.848 (°C) 

Form Figure 38 it can be seen that the post run calibration matches well with 
the data acquisition calibration curve. 
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5. Platinum Wire 5 

From Figure 39 it is clear to see the shift in calibration curves between the cal- 
ibration bath data and initial inplace calibration data. The calibration data has been 
altered while it was being inserted into the test chamber obviously. The initial data 
calibration was used for the experimental run calibration curve. The resulting cali- 
bration was: 



T s = 53.085 x R s - 260.049 (°C) 

From Figure 40 it can be further seen that a considerable error has been intro- 
duced into wire 5's calibration, an error of 3.7° C. This places an uncertainity value onto 
wire's 5 temperature data of at least 3.7° C, an unacceptable value. 
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Figure 30. Thermocouple Uncertainly: Uncertainly of thermocouple's 

tempertaure measurements as compared to a platinum resistance tem- 
perature probe. 
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Figure 31. Calibration of Platinum Wire 1: Calibration curve of wire 1 of Cali- 

bration Bath Temperature Increasing, Calibration Bath Temperature 
Decreasing, and of the Initial in place calibration data. 
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WIRE 1 CALIBRATION 




Figure 32. Calibration of Platinum Wire 1: Calibration curve of wire 1 of Initial 
in place calibration curve, and of Final in place calibration curve 
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TEMPERATURE (DEG C) 



WIRE 2 CALIBRATION 
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Figure 33. Calibration of Platinum Wire 2: Calibration curve of wire 2 of Cali- 
bration Bath Temperature Increasing, Calibration Bath Temperature 
Decreasing, and of the Initial in place calibration data. 
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Figure 34. Calibration of Platinum Wire 2: Calibration curve of wire 2 of Initial 

in place calibration curve, and of Final in place calibration curve 
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WIRE 3 C ALI BRATION 




Figure 35. Calibration of Platinum Wire 3: Calibration curve of wire 3 of Cali- 
bration Bath Temperature Increasing, Calibration Bath Temperature 
Decreasing, and of the Initial in place calibration data. 
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TEMPERATURE (DEG C) 



WIRE 3 CALIBRATION 
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Figure 36. Calibration of Platinum Wire 3: Calibration curve of wire 3 of Initial 
in place calibration curve, and of Final in place calibration curve 
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TEMPERATURE (DEG C) 
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Figure 37. Calibration of Platinum Wire 4: Calibration curve of wire 4 of Cali- 

bration Bath Temperature Increasing, Calibration Bath Temperature 
Decreasing, and of the Initial in place calibration data. 
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WIRE 4 CALIBRATION 



o 




Figure 38- Calibration of Platinum Wire 4: Calibration curve of wire 4 of Initial 

in place calibration curve, and of Final in place calibration curve 
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TEMPERATURE (DEG C) 
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Figure 39. Calibration of Platinum Wire 5: Calibration curve of wire 5 of Cali- 
bration Bath Temperature Increasing, Calibration Bath Temperature 
Decreasing, and of the Initial in place calibration data. 



62 



TEMPERATURE (DEG C) 



WIRE 5 CALIBRATION 




Figure 40. Calibration of Platinum Wire 5: Calibration curve of wire 5 of Initial 

in place calibration curve, and of Final in place calibration curve 



F. SUMMARY 

The following points were brought out during the calibration process. 

1. The process of moving the platinum wires from the calibration bath and installing 
them into the test chamber resulted in a shift in the calibration curves. 
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TEMPERATURE (DEG C) 



2. The most probable cause of all the calibration shifts were the lead wires in the back 
of the insert board being pulled and adding an additional tension onto the platinum 
wires, thereby stretching them. The fact that the calibration curves consistently 
shifted to an increase in resistance supports this theory. 

3. The maximum difference observed between the Platinum Wire Temperature Probe 
and the Type T thermocouples were 0.4° C. Therefore, a temperature uncertainity 
for the in place calibration will be set at 0.4° C. 

4. The differences between the wire's data acquisition calibration curve and the wire's 
fmal calibration curve did not exceed the thermocouples uncertainity in temper- 
ature. 
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APPENDIX B. SAMPLE CALCULATIONS 

In these sample calculations the output value of wire 1 will be used. 

A. DETERMINATION OF DERIVED VALUES 

1. Determination of the Average Bulk Temperature 

The average bulk temperature of the liquid w r as obtained by the arithmetic av- 
erage of the four immersed themocouple outputs. 

(T, + T 2 + T 3 + T 4 ) 

'■bulk 4_q 

(55.835 + 55.884 + 55.988 + 55.932) 



4.0 



= 55.9098°C 



2. Determination of Wire Current 

The current flow through the platinum wire is obtained by measuring the volt- 
age drop across a precision 2 ohm resistor. 



I 



Res 



PtWire ' 



^ 20 . 

(0.74605Volts) 

(1.99940hm) 
= 0.2381Ohms 



B. DETERMINATION OF PLATINUM WIRE RESISTANCE 

the resistance value of the platinum wires are obtained by measuring the voltage 
drop across the desired wire and dividing it by the current flow through it. 



^ PtWire ~ i 



PtWire 



PtWire 

(1.7233Volts) 
“ (0.2831 Amps) 
= 6.0874 Ohms 



1. Determination of Wire Surface Temperature 

The wire's surface temperature are obtained from the calibration curves shown 
in APPENDIX A, wire 1 will be shown as an example. 
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^ PtWlre 1 — 53.902 X R Ptwire ~ 254.470 
= 53.902 x 6.0874 - 254.470 
= 73.356°C 

2. Determination of Input Power 

The input power of the wires are determined from the current flow through the 
wire times the voltage drop across the wire. 

Power] = pnvire\ ^ Pnvire 1 

= 1.7233Volts x 0.283 lOhms 
= 0.48787Watts 

3. Determination of Heat Flux 

The heat flux dissipated by the platinum wire is determined from the inputted 
power divided by the wire's surface area. 



HfluX] = 



Power, 

Ai 

(0.48787Watts) 

1.5585E-5m 2 



= 31,303.44Watts/m 2 



where, 



A, = n x D, x L, 

= n x 5E-5 m x 0.09922m 
= 1.5585E-5'm 2 

C. DETERMINATION OF FLUID PROPERTIES 
1. Film Temperature 

The film temperature is the arithmetic average of bulk fluid temperature, and the 
wire's surface temperature. 



_ bulk+T PtlVire) 

'film ~ 2.0 

(54.0°C + 58.0°C) 

2.0 

= 56.0°C 
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2. Thermal Conductivity 

From Figure 5 of the 3M Corporation Flourinert Product Manual [Ref. 3 : p. 
15], the thermal conductivity has been determined to be: 

(0.6033 -0.001 15 xTyj /m ) 

k= w 

= 0.05389 Watts/m°C 



at T„ /m = 56°C 

3. Liquid Density 

Using the expressions in Table 4B and the constants presented in Table 4C of 
the 3M Corporation Product Manual, [Ref 3: p.10], the liquid density has been deter- 
mined to be: 



P/ = (1.740 -0.00261 x Tyi /m ) x 1000 
= l,593.8kg/m 3 



at ^ film — 56 C 

4. Kinematic Viscosity 

Using Figure 3 in the 3M Corporation Product Manual, [Ref 3: p.13] and de- 
termining an expodintial curve fit valid from 0°C to 90°C yields: 

v = A exp( t ^ + 273 . 15 x ) 

A = 1.203952E — 8m 2 /s 
B= 1,058.4109°K 



at T fum — 56 C 



v = 3.000E - 7m 2 /s 



5. Specific Heat 

From Figure 4 of the 3M Corporation Flourinert Product Manual, [Ref. 3 : 
p. 14], the specific heat has been determined to be: 

C P = (0.241 1 1 1 + 3.70374E - 4 x T fllm ) x 4186 
= 1,096. lJ/kg°C 



at Tfl /m = 56° C 
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6. Thermal Expansion Coefficient 

Using the expression in Table 4B and the constants presented in Table 4C of the 
3M Corporation Product Manual, [Ref. 3: p. 10], the thermal expansion coefficient has 
been determined to be: 



» 0.00261 

P (1.740 —0.00261 x Tfii m ) 
= 0.00164 1/°C 



at T„, m = 56°C 

7. Thermal Diffusivity 

Is determined to be at = 56°C 



p x C p 

= 3.0848E - 8m 2 /s 



8. Prandlt Number 

Prandlt Number is determined to be at T^ = 56°C 



Pr = 



_v 

a 



= 9.725 
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APPENDIX C. UNCERTAINTY ANALYSIS 



A. UNCERTAINTY IN SURFACE AREA 
co D = 0.003 mm 

o) L = 0.8 mm 

A = n D L 
where: 

D = Wire Diameter 
L = Length of platinum wire 
A = Surface Area of platinum wire 

dA , T s 

Id " ( ” L) 

dA , p.x 

1l ' ( ” D) 

Therefore, 



co A = 


1 1 

IT 

r* 

S 

to 


+ (n 


D coj 2 ] 1 ' 1 






( Ud > 


M 


f Vll 


A 




A D , 


V L ) _ 




[( 


' 0.003 > 


M 


' 0.8 Y 




l 


, 0.05 ) 


v 99.22 ) 



= 0.061 or 6.1% 

A 

c o A = 9.43529E -7 m 1 
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B. UNCERTAINTY IN POWER 



co r2ii = 0.008% 4- 8/xV [Ref 13, /?g.587] 
cov.pimre = 0.008% + 8/cV [Ref 13,/>g.587] 



(o R2Sl = ± 0.00 Id max value 



Q — ^ PtWire x I 



P t Wire 



= V 



'2Q 



Pt Wire 



2 Q. 



dQ 



V 



2a 



dV 



Pi Wire 

dQ 



R. 



7Q 



' Pt Wire 



av. 



2 a 



R 



2 a 



JQ_ 

dR 



= -V 



v. 



2 a 



2 a 



PtWire -r 2 

K 2a 



<y<2 = 



V 2a .. \l . ( V p t wir e 

£ X 03 V, PtWire j + \ 



X CO 



v 2a 



V,2£2 



+ -V 



PtWire x 



2 a 



X CO 



2 a 



2 a 



1/2 



CO/ 



MV, PtWire \2 ( ^V.lQ. \ 2 



'PtWire 



+ 



+ 



' 2 a 



t0 2a 
R. 



2 a 



1/2 



The worst case of uncertainty will occur when the voltages measured for V PtW , r , and V 2Q 
are small, less than 0.3 Volts. 



CO/ 






(0.008% + 0.0027%) 2 + (0.008% + 0.0027%) 2 + 



(t 



0.001 0/o \2 



9994 x 100 



O’ 



1/2 



= 0.0524% or 5.2384E - 4 
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C. UNCERTAINTY IN HEAT FLUX 



dq _ 1 

dQ A 

dq_ = 

dA a 2 




= [(0.06 1) 2 + (5.2384E-4) 2 ] 1 ' 2 



= 6.1002E-2 or % 6.1 

D. UNCERTAINTY IN TEMPERATURE 

The uncertainty in the thermocouple measurements (<y rc ) is determined by the un- 
certainty in the platinum wire temperature probe (c o PtPrQb , ) and the error in the 
thermocouple's in matching the probe temperature (co Arc ). 

W TC = [(^PtProbe ) 2 + ( W A TC ?]' 12 

= 0.400 1°C 

E. UNCERTAINTY OF WIRE SURFACE TEMPERATURE 
From the calibration formula 



co 



0.1 



Q. 



c o To = 0.4001 °C 
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Trar/ “ £ Rprmre + T 0 



R 



20 



~ PtWire y + T 0 + <5T 



20 . 



Where: 



e = slope of calibration curve 



T 0 = axis intercept for R = 0 



V ptw, re = Voltage drop across Platinum Wire 



R 2Q = Resistance of in series Resistor 



W 2 o = Voltage drop across series Resistor 

<5T = Difference between Data Acquisition and Final Calibration for a particular data po 



3 1 



surf 



v 



Pt Wire 



3l 



R, 



v ~ 

v 20 



31 



surf 



5Tn 



= 1 



5T 



surf 



351 



= 1 



31. 



surf 



R 



av 



= E 



2a 



Ptmre 



f 2Q 



dl 



surf 



5R 



= E ' 



1 PtWire 



2 a 



'2a 



dl, 



surf 



5V 



= -cV 



R 



20 



20 



PtWire 2 

^20 



12 



0J Tsurf 




V Pi Wire 



( 




/ + 0 "ar) 2 1/2 



Using data from Wire 1 at the following point, 
W PtWir , = 0.7 1673 Volts 
V 2n = 0.24529 Volts 
R 2Q = 1.9994 Ohms for wire 1 

«.*.« = 53.902°C /Q 
<5T = 0.02°C 

This results in an uncertainty in temperature of, 

w Tsurf.wire 1 = 0 . 8 54 ° C 
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APPENDIX D. DATA ACQUISITION PROGRAMS USED 



A. PLATINUM WIRE RESISTANCE CALIBRATION PROGRAM 

The following program was used to obtain the Resistance value's for the platinum 
wire specimans used in chapter 2. This program is written in HP basic language for use 
with controlling the HP-3852A Data Acquisition unit. 



10 DIM Ohms (0:61 ) 

20 OUTPUT 709 ; " INTE6ER I" 
30 OUTPUT 709 i " RST 000" 

40 OUTPUT 709 ; "US^, 000" 

50 OUTPUT 709 ; "CONF OHMF “ 
60 OUTPUT 709; "NPLC 1 .0" 
70 OUTPUT 709; "NRDGS 62" 
80 OUTPUT 709; "DELAY 0.1" 
90 OUTPUT 709; "RANGE 10" 



100 


OUTPUT 709; "I = 


99" 


1 1 0 


OUTPUT 709 ;"I = 


1 +1" 


120 


OUTPUT 709 ; " MEAS OHMF . I" 


130 


ENTER 709;Ohms( * ) 


140 


A 1 o= 1 . E+30 




150 


J=1+J 




160 


Ah l = 1 .E-30 




170 


A=0 . 




1 80 


FOR K= 1 TO 60 




1 90 


A=A+Ohms( K ) 




200 


IF Ohms ( K XAl o 


THEN 


210 


Alo=Ohms ( K ) 




220 


END IF 




230 


IF Ohms ( K ) > A h l 


THEN 


240 


Ahi=Ohms( K ) 




250 


END IF 




260 


NEXT K 




270 


B=A/60 . 




280 


PRINT "WIRE SAMPLE ";J 


290 


PRINT "60 RESISTANCE AVERAGE IS “;B 


300 


PRINT "10 RESISTANCE MEASURED IS ";Alo 


310 


PRINT "HI RESISTANCE MEASURED IS ";Ahi 


320 


IF J=5 THEN 




330 


GOTO 360 




340 


END IF 




350 


GOTO 1 10 




360 


GOTO 370 




370 


DIM T empsl 0:61) 




380 


K = 0 . 




390 


OUTPUT 709; "RST 


000" 


400 


OUTPUT 709; "USE 


000" 


410 


OUTPUT 709; "INTEGER I" 
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420 OUTPUT 709 ; "CONF TEMPT" 

430 OUTPUT 709; "NPLC 1.0” 

440 OUTPUT 709; "NRDGS G2" 

450 OUTPUT 709; "DELAY 0.1" 

460 OUTPUT 709 ; "RANGE AUTO" 

470 OUTPUT 709; "1=299" 

480 OUTPUT 709; "1=1+1" 

490 OUTPUT 709 ; " MEA5 TEMPT, I" 

500 ENTER 709 ; Temps ( » ) 

510 K=K+ 1 
520 A=0. 

530 FOR L= 1 TO 60 
540 A=A+Tenps(L) 

550 NEXT L 
560 Terop=A/60 . 

570 PRINT "THERMOCOUPLE" ;K; "TEMPERATURE" ; Temp 

580 IF K=4 THEN 

590 GOTO 620 

600 END IF 

610 GOTO 480 

620 END 



75 



B. MAIN DATA ACQUISITION PROGRAM 

The following program was used to obtain the Temperature value's for the platinum 
wire specimans used in chapter 3. This program is written in HP basic language for use 
with controlling the HP-3852A Data Acquisition unit. 



10 
20 
30 
40 
50 
G0 
70 
80 
90 
100 
1 10 
120 
130 
140 
150 
160 
170 
1 80 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 



DIM Temps ( 0: 16 ) ,Tk< 1 : 7 ) ,Tt ( 1 : 1 00 ) ,Flux( 1 : 1 00 ) ,Tb ( 1 00 ) 

L=0 

J=0 

L=l+ 1 

IF L=8 1 THEN 
GOTO 1450 
END IF 

OUTPUT 709r ,v RST 000" 

OUTPUT 709; "USE 000" 

OUTPUT 709;" INTEGER I" 

OUTPUT 709 ; "CONF TEMPT" 

OUTPUT 709 ; "NPLC 1.0" 

OUTPUT 709 ; " NRDGS 17- 
OUTPUT 709; "DELAY 0.1" 

OUTPUT 709; "RANGE AUTO- 
OUTPUT 709; "1=299" 

OUTPUT 709; "1 = 1 + 1 ” 

OUTPUT 709 ; " MEAS TEMPT, I" 

ENTER 709 ; T emps ( * ) 

J = J+1 
A=0 . 

FOR K= 1 TO 15 
A=A+Temps( K ) 

NEXT K 
Tk< J ) = A/15. 

IF J=7 THEN 
GOTO 300 
END IF 
GOTO 170 

DIM Uol t s ( 0 : 1 1 ) ,Uk( 1:10) 

J=0. 

Ct=0. 

OUTPUT 709 ; "RST 000- 
OUTPUT 709; "USE 000" 

OUTPUT 709; "INTEGER I" 

OUTPUT 709; "CONF DCV" 

OUTPUT 709; "NPLC 1.0" 

OUTPUT 709; "NRDGS 12" 

OUTPUT 709; "DELAY 0.1" 

OUTPUT 709; "RANGE AUTO" 

OUTPUT 709; "1 = 99" 

OUTPUT 709; "1 = 1 + 1 “ 

OUTPUT 709; "MEAS DCV, I" 

ENTER 709 ; Volts( * ) 

J = J+1 
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4G0 

470 

480 

480 

500 

510 

520 

530 

540 

550 

550 

570 

580 

590 

600 

510 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 



IF J=5 THEN 
GOTO 420 
END IF 
Ct=Ct+l 
A=0 . 

FOR K= 1 TO 10 
A=A+Ool ts< K ) 

NEXT K 

Ok < Ct >=A/ 1 0 . 

IF J=6 THEN 

GOTO 590 

END IF 

GOTO 420 

DIN Ooltw( 0: 1 1 ) 

J=0. 

Ct=5 . 

OUTPUT 709?' RST 000" 

OUTPUT 709; "USE 000" 

OUTPUT 709; "INTEGER I" 

OUTPUT 709 ; "CONF DCO" 

OUTPUT 709; "NPLC I .0" 

OUTPUT 709 ; "NRDGS 12" 

OUTPUT 709; "DELAY 0.1" 

OUTPUT 709; "RANGE AUTO- 
OUTPUT 709; "1=109" 

OUTPUT 709; "1=1 + 1 " 

OUTPUT 709; "NEAS DCO , I " 

ENTER 709 ; Ool tu( * ) 

J = J+1 

IF J=5 THEN 
GOTO 710 
END IF 
Ct=Ct+l 
A=0 . 

FOR K= 1 TO 10 
A=A+Oo 1 tw( K ) 

NEXT K 

Ok < Ct > = A / 1 0 . 

IF J=6 THEN 
GOTO 880 
END IF 
GOTO 710 

Bt=( Tk< I >+Tk(2)+Tk(3)+Tk<4>>/4. 
01 =0k. ( 6 ) 

1 1 =0k ( 1 )/l .9994 
02=0k( 7 ) 

I2=0k( 2 )/2 .0076 
03=0k< 8 ) 

I3=0k< 3 )/ 1 .9997 
,04=0k( 9 ) 

1 4=0k ( 4 )/2 .001 8 
05=0k ( 10 ) 

I5=0k( 5 )/2 . 0025 
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990 P1=I1*U1 

1000 P2=I2*U2 

1010 P3=I3*U3 

1020 P4=I4*U4 

1030 P5=I5*U5 

1040 R 1 =U 1 / 1 1 

1050 R2=U2/I2 

1 0B0 !R3=U3/I3 

1070 ! R4=U4/ 1 4 

1080 !R5=U5/I5 

1090 Flux 1 =P 1 / 1 . 5585284E-5 

1100 Flux2=P2/l .5585284E-5 

1110 Flux3=P3/1 .5710005E-5 

1120 Flux4=P4/1 .54G05G3E-5 

1130 Flux5=P5/1 .5710005E-5 

1140 T1 =( 53 . 902 1 23*R 1 )-254. 46995 

1150 T2=< 53 . 666456*R2 )-253. 00959 

1150 T3=( 53 . 1 4 1 047*R3 )-252 . 7291 3 1 

1 170 T4=< 53. 483646*R4)-252. 848797 

1180 T5=(53.085508*R5)-260. 04883 

1190 Tc = ( Tk< 5 )+Tk( 6 ) )/2 . 

1200 PRINT ’ TEMPERATURES" 

1210' PRINT " AMBIENT AIR BULK FLUID CONDENSOR SURFACE" 

1220 PRINT Tk ( 7 ) ,Bt ,Tc 
1230 PRINT “ 

1240 PRINT “WIRE . TEMPERATURE CURRENT UOLTS * 

1250 ! PRINT M (AMPS) ( UOLTS ) (DEGREE C) 

1260 PRINT “ 1 " , T 1 ,11 ,V1 

1270 PRINT "2” ,T2 ,12 ,V2 

1280 IPRINT "3“ ,T3 , 13 ,U3 

1290 ! PRINT "4" ,T4 ,14 ,U4 

1300 IPRINT "5“ ,T5 ,15 ,U5 

1310 PRINT “ HEAT FLUX POWER RESISTANCE 

1320 IPRINT " (WATTS) (OHMS) ( WATTS/M2 ) 

1330 PRINT “ 1 “ .Flux 1 ,P 1 ,R 1 

1340 PRINT “2" ,Flux2 ,P2 ,R2 

1350 IPRINT "3" ,FLUX3 ,P3 ,R3 

1360 IPRINT "4” ,Flux4 ,P4 ,R4 

1370 IPRINT “5" ,Flux5 ,P5 ,R5 

1380 PRINT " " 

1390 Tt ( L )=T 1 
1400 Tb(L )=Bt 
1410 Flux( L )=Flux1 
1420 WAIT 100 
1430 GOTO 30 
1440 GOTO 1530 
1450 M= 1 

1 4G0 PRINT " WIRE TEMP HEAT FLUX BULK TEMP" 

1470 IF M=3 THEN 
1480 GOTO 1530 
1490 END IF 

1500 PRINT M ,Tt ( M ) ,Flux( M > ,Tb( M ) 

1510 M=M+1 

1520 GOTO 1470 
1530 END 
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C. NUSSELT EVALUATION PROGRAM 

The following program was used to obtain the Nusselt Number for the platinum 
wire specimans used in chapter 3. This program is written in Fortran language for use 
with controlling the IBM-3803 Main Frame Computer. 



FILE: NUSSELT FORTRAN A1 



SUBROUTINE NUSSELT (X1,T1,M,Q) 

C 

C THIS SUBROUTINE WILL CALCULATE THE HEAT TRANSFER COEFFICIENT FOR THE 
C DATA OBTAINED FROM THE EXPERIMENTAL RUNS BY USING THE CHURCHILL AND 
C CHU CORRELATION 
C 

DIMENSION X1(X),T1(X),Q(X) 

REALX4 X1,T1,Q,K,NU,NUD 

OPEN ( 29 , FI L E= 1 RAL Y ' , STATUS= * UNKNOWN 1 , F0RM= * FORMATTED • ) 

C 

C DIAMETER OF THE WIRE BEING ANALYSED IS NOMINALLY .05 MM 
C 

D=5.0E-5 

C 

C GRAVITY TERM OF 9.81 M/S2 
C 

G=9 . 81 
C 

DO 100 1=1, M 
C 

C FILM TEMPERATURE OF THE PHYSICAL PROPERTIES 
C 

CCC TF = (XI ( I ) + Tl(I) + 56 . 0 )/2 . 

TF = (XI ( I ) + 2.xTl(I))/2. 

C 

C TEMPERATURE DIFFERENCE BETWEEN THE WIRE AND BULK FLUID 
C 

CCC TD =(X1(I) - Tl(I) + 56.0) 

TD = XI (I) 

C 

C THERMAL CONDUCTIVITY OF FLUID W/MK 
C 

K = (0.6033 - 0.00115*TF)/10. 

C 

C BULK MODULUS OF THE FLUID 1/K 
C 

BETA = 0.00261/(1.740 - 0.00261XTF) 

C 

C KINEMATIC VISCOSITY OF THE FLUID M2/S 
C 

NU = 1.203952E-8 x EXP( 1 058 . 41 09/ ( TF+27 3 . 15 ) ) 

C 

C SPECIFIC HEAT OF FLUID 
C 

CP = (0.24111 + 3.70374E-4 x TF)*4186.00 
C 

C FLUID DENSITY KG/M3 
C 

RHO = (1.740 - 0.00261*TF)xl000.0 
C 
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oooo ooo ooo ooo on 



THEMAL DIFFUSIVITY 
ALFA = K/(RH0*CP) 

PRANDLT NUMBER OF FLUID 
PR = NU/ALFA 

RALYEIGH NUMBER OF FLUID MUST BE 10**-5< RA < 10**12 

RA = (G*BETA*TD*D**3 . )/(NU*ALFA) 

PRINT *, RA 

NUSSELT NUMBER OF FLUID NU=H*D/K 

TERM = ((0.518*CRA**C0.25))*(l+(0.559/PR)**0.6)**C-5./12.))X*15 
|+(0.1*RA**(l./3.))x*15.)**(l./15.) 

NUD=2.0/CL0G(1.+ 2 ./TERM) ) 

I 

THERORICAL HEAT FLUX WATTS/M2 



QCI) = K*NUD*TD/D 

c 100 WRITEC29,500) TD, RA, NUD, PR, Q(I ) 

500 FORMAT ( 5( F12 . A, 2X) ) 

CLOSE (29) 

RETURN 

END 
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